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We propose a new scheme for parallel spatially multimode quantum memory for light. The scheme 
is based on counter-propagating quantum signal wave and strong classical reference wave, like in a 
classical volume hologram, and therefore can be called a quantum voiume hologram. The medium 
for the hologram consists of a spatially extended ensemble of atoms placed in a magnetic field. The 
write-in and read-out of this quantum hologram is as simple as that of its classical counterpart 
and consists of a single pass illumination. In addition we show that the present scheme for a 
quantum hologram is less sensitive to diffraction and therefore is capable of achieving higher density 
, of storage of spatial modes as compared to previous proposals. A quantum hologram capable of 

^ ■ storing entangled images can become an important ingredient in quantum information processing 

' and quantum imaging. 
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PACS numbers: 42.50.Ex, 03.67.Mn, 37.10.Jk, 42.40.-i 

I. INTRODUCTION 



i-C ' Quantum memory for light is an essential part of quantum information protocols, such as quantum repeaters, 
^H, distributed quantum computation, and quantum networks. A number of approaches based on storage in atomic 
ensembles were developed recently, including quantum non-demolition (QND) interaction, electromagnetically induced 
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transparency (EIT), Raman interaction and photon echo. For a comprehensive recent review on quantum interfaces 
between light and matter see Multimode quantum memories are in the center of current research due to their 
potential for enhanced storage capacity and state purification and error correction, e.g., in quantum repeaters (2j. 
— 1 . Spatially-multimode parallel quantum protocols for light without memory have been central for the field of quantum 
imaging. Examples are quantum holographic teleportation d, 0] and telecloning Q , and quantum dense coding of 
optical images 6j . The spatially multimode light in an entangled Einstein - Podolsky - Rosen (EPR) quantum state 
Q has been recently experimentally demonstrated via four- wave mixing in Q. 
£S| \ A quantum memory protocol based on a QND interaction, for the single-mode case experimentally demonstrated 
I/"") ■ in 0, achieves storage in the ground state of an ensemble of spin-polarized atoms. An extension of the QND-based 
scheme of quantum memory to spatially multimode configuration, able to store quantum images, has been proposed 
recently in '10']. In this proposal two passes of light are required to store, and another two - to retrieve quantum state 
of an image from the hologram. In addition, perfect squeezing of the initial state of light and/or atoms is required 
ON . for an ideal performance of the holographic memory. 

Another approach to the multimode quantum memory based on phase-matched backward propagation retrieval 
out of the EIT and the Raman type memories has been recently proposed for the multiple frequency-encoded qubits 
storage in a single ensemble . A gradient echo memory has been proposed for storage of several frequency encoded 
modes fl2| . 

The optical image storage has been demonstrated in recent experiments [H, [3 [l5| . There was observed [r| a 
several nanoseconds delay of optical pulses which contain on average less than one photon and carry two-dimensional 
images. The delay has been achieved by using slow light in an atomic ensemble. However, preservation of quantum 
field properties was not demonstrated in the paper. The storage of classical images in a warm atomic vapor by EIT 
interaction has been investigated in [3, EH • 

In this paper we present a new scheme of the multimode quantum memory which we call a quantum volume 
hologram. Our proposal makes use of two concepts. The first one stems from a volume hologram, proposed by 
Denisyuk [Tfjj for classical storage of optical images. The volume hologram is written by the counter-propagating 
signal and reference waves. There are two sublattices, produced by the waves interfering in the medium, each of them 
stores one quadrature of the signal field. Since both quadratures are stored, there are no virtual and real images 
during the readout, by contrast to a thin hologram. 

The second concept comes from [l7j where it was shown that a combination of a constant magnetic field with the 
QND interaction allows to couple two components of the collective atomic spin used for the quantum memory to two 
quadratures of light . 

In the present paper we propose to write volume hologram (which implies the counter-propagating geometry) onto 
an atomic ensemble with spins rotating in a constant magnetic field. Hence, in our model all degrees of freedom of 
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FIG. 1: Schematic of the level structure and of the write (a) and readout (b) stages of quantum volume hologram 



atoms and light are symmetrically involved in the interaction. We show that an input state of spatially-multimode 
light can be written onto the quantum volume hologram in a single pass. A sideband (with respect to the carrier 
frequency of the strong field to) +f2 spectral component of light is written onto the collective spin coherence wave, 
which propagates in the medium with a certain phase velocity. 

Our analysis shows, that the volume hologram for the ground state spin J = 1/2, when considered in the limit of 
a single spatial mode, has features similar to those described in [l8| for A-scheme of atomic levels (possible only for 
J > 1/2). It appears that an effective increase in the number of the system degrees of freedom due to the counter- 
propagating geometry in combination with the spin rotation has an effect similar to the increased complexity of the 
atomic level structure in the A-scheme based memories. 

The paper is organized as follows. First, we discuss the basics of the model and derive the light and matter 
equations of motion in the paraxial approximation using averaging over fast oscillations in space and time. Next, 
we consider transfer of an input multimode quantum state encoded in the signal wavefront at the write and readout 
stages. Finally, the stored transverse modes number is estimated for the volume and the thin [10j quantum hologram 
and we come to conclusion. 



II. SINGLE PASS VOLUME HOLOGRAM WITH SPATIAL RESOLUTION 

The scheme, illustrating our quantum memory protocol is shown in Fig. [TJ We consider an ensemble of motionless 
atoms which we for simplicity assume to have an angular momentum of J = 1/2 both in the ground and in the excited 
states, located at random positions. The long-lived ground state spin of an atom J a is initially oriented along the 
constant magnetic field in the vertical direction x. The atomic spins rotate around the vertical axis with a circular 
frequency Q. A classical off-resonant rr-polarized plane wave at frequency luq with a slowly varying amplitude A x 
(assumed to be real) propagates in —z direction. The input signal wave is a weak quantized y-polarized field at the 
same frequency u>o, propagating in +z direction. In what follows we consider this multimode input field with a slowly 
varying amplitude A y (r, t) <C A x in the paraxial approximation. 

In order to construct the Hamiltonian of the system we first consider an atomic layer with the thickness (in z 
direction) smaller than A. Within this layer the phase difference between the signal field and the driving wave is 
constant. Hence, the interaction within the slice is described by the well known quantum non-demolition (QND) 
Hamiltonian [l[. The QND light-matter interaction in the each layer leads to two basic effects: (i) the Faraday 
rotation of light polarization due to longitudinal quantum z-component of collective atomic spin of the slice; and (ii) 
the atomic spin rotation, caused by the unequal light shifts of the ground state sub- levels with m z = ±1/2 in the 
presence of quantum fluctuations of circular light polarizations. The relevant part of the Hamiltonian is [l[: 



J dr^ J2®S,(f, t)S(r- r a ). (1) 



UJ e g — UlQ Jy 

Here aj eg is the frequency of the atomic transition, d is the dipole matrix element, and fco = ujqc. For the counter- 
propagating signal and driving waves, the ^-component of the Stokes vector S z (r,t) = 2A X Im \A y (r, t)e 2k ° z ~\ in ((T|) 
rapidly (i.e. on the scale A) oscillates along z. The amplitude A y is defined via 

dk z f dq 



A y (z,p,t)= J — ^ J j^-^ y/oj(k)/k a y (k)exp[i(q- p+ (k z - k Q )z - (u(k) - w )t)], 
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here a v (k) and a y (k) are the annihilation and creation operators for the wave k, which obey standard commutation 

relations [a y (k), al(k')] = {2ir) 3 8{k — k'), [a y (k),a y (k')] = 0. By using these commutation relations in the paraxial 
approximation, one finds [7|] the commutation relation for the slowly varying amplitude of the quantized signal field 

Ay, 

[Ay(z, p, i),4(z', ?',()] = cS(f-f ) (2) 



where r — (p,z), k — (q, k z ). We do not consider here the y-polarized quantized field co-propagating with the driving 
wave because its evolution is independent of the signal wave under consideration. 

We introduce the density of the collective spin as J(r) = ^ a J a 8(r — r a ). The averaged over random positions of 
the atoms commutation relation for the y, z components of the collective spin is 

[J y (r),J z (r')} =iY J {J a x )5{r-r a )5{r'-r a ) a = in a {J a x )5{r- ?). 

a 

Here n a is the average density of atoms. The field-like variable for the spin subsystem, 

B{r,t) = {J y (r, t) + U z (r, t)}/ V2n„< J«), 
obeys the standard boson commutation relation: 

[B(r,t),B^r',t)]=S(r-r'). (3) 

In view of the fact that the atoms are prepared in the spin up state, the operator B^ can be considered as the 
creation operator for atoms in the spin down state, or the collective projection operator | + 1/2) — > | — 1/2). 

The full Hamiltonian of our model includes the energy of a free electromagnetic field, the interaction of the atomic 
spin with constant magnetic field and the effective Hamiltonian of the QND interaction. The Hamiltonian reads, 

H = J df | ^ A\ {z, p, t) A y {z, p, t) + hi! fit ( Zj t) B{z, p, t)- (4) 



V2TTk Q \d\ 2 



tMJS) A x (z, t) [B(z, p, t) - h.c] [A y {z, p, ty 2kaz - h. 



We describe the evolution of our system in the Heisenberg picture. Using the commutation relations @ and ^ for 
the field and atomic variables, after simple transformations we obtain: 



d i „, 19,, 

Oz-W^^cm^^^ 



'LT 



Im_B(z, p, t)e 



-i2ko z 



(5) 



^B(z, p, t) = -itlB(z, p, t) + -J=Im [A y (z, p, t)e ak ° z ] . (6) 



Here T is the duration of the flat-top pulse of the driving field and L is the atomic cell length. The dimensionless 
coupling constant 



2irk \d{' 



k= U| 1 y/2n a {J°)A*LT, (7) 

h(UJ e g - IVO) 

should be of the order unity for the memory to work. The coupling constant can be written as k 2 = agij, where eta 
is the resonant optical depth and r\ is the probability of spontaneous emission [l[ . Since r\ <C 1 is required in order to 
neglect the effect of spontaneous emission from the excited atomic level, the usual condition for an efficient quantum 
interface ao — \ 2 n a L/2-K ^> 1 should be fulfilled. We introduce the Fourier transform via 

a(z, q, t) = [ dpA y (z, p, t)e~^ p ~, (8) 
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and similar for atomic variables, and arrive at the set of basic equations in the Fourier domain, 

' 2 1 d \ , _ k 1 



_5 

0* 



; jT 

2/cq c i9i 



a(z,q,t) 



\b{z,q,t) - b ] {z,-q,i)] e" 



i'2koz 



(9) 



— b(z,q,t) = -i£lb(z,q,t) 



<IJT 2i 



i[fl(z,Ci) e *^ fl t( z ,^() e - 2l » 2 ] 



(10) 



The considered above field and atomic amplitudes are rapidly oscillating at frequency of the order O. In what follows 
we introduce the slowly varying amplitudes of collective spin b±k(z, q, t) in the rotating at Q frame, and slowly varying 
amplitudes a±n(z, q, t) for the signal field sidebands at frequencies loq ± f2. 

Next, we should take into account that the field amplitudes are defined as slowly varying along the z axis, but 
the spin amplitudes are not, as seen from (|10p. The fast modulation of the collective spin at the longitudinal spatial 
frequency 2fc is just a consequence of the counter propagating geometry of the volume hologram. The thin atomic 
slices discussed above have a length of the order of a fraction of A. This imposes limitations on the atomic motion 
during the storage time: the atoms should not transport coherence to another slice. A solid state, an optical lattice, 
or an ultra cold atomic ensemble should be used to fulfil this condition. The slow in space and time variables look 
like 



a±u(z,q,t) = a(z,q,t)e 



±iOt 



(11) 



b ±k (z,q,t)^b(z,q,t)e^ 2k ° z \ 

Quantum memory is realized via the dynamical coupling between an(z, q, t) and bk(z, q, t). We insert these variables 
into $9§, (fT0|) and perform averaging on a time scale 1/S]<t m «T and on a spatial scale in z direction A <§C z av <C L. 
This implies that the pulse duration and the cell length are taken large enough, HT 3> 1 and koL 3> 1, where the last 
condition is typical also for classical volume holograms. 
The equations of motion read, 



—a n (z,q,t) = 



cdt \c 2ko 



a n (z,q,t) - i 



2VLT 



bk(z,q,t), 



(12) 



d_ 

dt 



bk(z,q,t) 



2VLT 



a n (z,q,t). 



(13) 



In what follows we neglect the retardation effects. Both the pulse length in space, estimated as c/5£l, where <5S1 is the 
signal spectral width, and the spatial length c/f2 of modulation at frequency Q are assumed to be much larger than 
the cell length, hence 5J7/c<C f2/c<C This allows to omit the terms ~ 1/c on the right side of (fT2"|) as compared 
to the term ~ 6/ s/LT, which is estimated as ~ a/L. Here we make use of ([T3|) and find for typical values k ~ 1 and 
t ~ T an estimate b ~ yjT/L a. 

Next, we go over to the amplitudes 

/ „2 \ / „2 \ 

(14) 



an(z,q,t) = an(z,q,t)exp i-^r 2 



b k (z,q,t) = b k (z,q,t)exp ( i-^r z 



and arrive at the equations 



— an{z,q,t) 
oz 



2VLT 



bk(z,q,t), 



Ft b k (z,q,t) 



2VLT 



an(z,q,t). 



(15) 



As compared to the equations, previously discussed (e.g., [3, [l9| and references therein) for the A-scheme Raman- 
type models of spatially single mode quantum memory, we do not encounter in (|15p the deleterious effect of the 
Stark shift of atomic levels by a classical control pulse. The level shift makes the phase matching and evolution 
more complicated, especially for a time- dependent profile of the control field. This simplification is due to the more 
symmetrical nature of the QND interaction. 

It follows from (|15|) that in our parallel quantum memory, based on the volume hologram, the signal and the spin 
coherence waves with an arbitrary q (in paraxial approximation) are interacting like the waves propagating along z 
direction in a spatially single mode Raman-type memory. The diffraction does not modify the state exchange between 
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light and matter, and our memory is able to store as many orthogonal spatial mode with different q as one can fit into 
a hologram of transverse area S. We take Aq x y ~ 2ir/yS, \q\ <C 27r/A, and assume that the number of transverse 
modes is much less than 4S/X 2 . 

The solution in terms of the amplitudes (|14[) is found as an extension to the spatially multimode case of the results 
previously obtained for the Raman-type model, 



a Q (z,q,t) = a n (0,q,t)--^= [ dt J h I ) 5n(0, q, t') (16) 



2VLT Jo V t-V IV LT 



—i- 



2VLT Jo \ V LT 



dz' J M — — — b k (z',q,0), 



b k (z, q, t) = b k (z, q, 0) - — J= J dz'^j j^— Ji f «y (Z £ / ^ I M* 7 , <f, 0) (17) 



-£= / dt' Join 
/If Jo \ 



Consider the dimensionless coordinates and amplitudes, 

^ = 7", an(z,q,t)VT= a n (£,q, r), b k (z,q,t)VL = /3 k (^,q,T). (18) 

Here a(z, p, i)^a(z, p, i) corresponds to the number of signal photons per cm 2 of the beam cross section during the 
interaction time T, and /3(z, /?, t)^/3(z, p, i) gives the number of flipped spins per cm 2 of the hologram cross section at 
the length L. The equations above become 

4T%»= f dr'G l {T-T',l)a < t\q,T , )-i f d£G {l r)fif l) (£', q), (19) 
Jo Jo 

p k out) ^q) = [ dt'Giit-eAM^ie^-i [ dT'G (i-T',o& ( n n) (qy), (20) 



where ag n \q,r) = a n (0,q,r), ag ut \q, T ) = a a (l,$,T), and ^(tq) = &(£,£0), P {out) {^q) = The 
integral kernels are given by 

Go(p,q) = 2 J o(«Vw)' Gi(p, g) = J(p) - -Ji(Ky/qp)^/q/pd(p), (21) 

with the nth Bessel function of the first kind denoted by J n . 

Similar to the analysis of [TH, we make use of the fact that the kernels Go,i share eigen functions, 

Go0—x,y) = ^MyjXiUl-x), (22) 

i 

Gi(y-x,l) = ~^(j>i{y)ni<j>i(l-x), 
and their eigen values satisfy the constraint Af + /if = 1 for all i. By applying the variables decomposition of the form 

4 m) (<?» = Eag^aWiCi-T), (23) 

3 
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FIG. 2: The first three eigen functions: cj>\, 02, and <f>3 (the bold, thin and dashed line respectively) for k = 4. The eigen values 
are Ai « 0.988, A 2 » -0.518, and A 3 » 0.043. 



one arrives at the input-output transformation for one cycle of the light-matter interaction: 



~ {out) 



(24) 



The transformation is of a beamsplitter type, it is unitary and hence respects the correct bosonic commutation 
relations. The Eq. (j24|) shows, that if the light signal stored in an atomic hologram has the temporal profile of an 
eigen mode of Go with the eigen value close to unity, the initial state of the atomic spins in the mode i in erased 
since the corresponding eigen value fii is close to zero. 

In order to read out the image stored in the hologram, one has to repeat the procedure and pass an- 
other classical light pulse through the atoms. To describe the whole write-readout cycle of parallel quantum 



memory, we combine the transformations 



{a$ in \q,T)J* in \t,q) 
(see (HI])). This yields, 



{ 



' ' '* } (<f> T ) = ex P ( ~ i ^ L J \ Mi a n, ( r ) 



~W{in) 
-R(out) 



(q, r), /3^ m ^(£, q)\ — > A^ (otlt) (£, (?) for the write stage, and 
(q, t) for the readout stage, and restore the diffraction factor 



{q)-i\iYj., 



nWUn) , -* 



iX r 



W(in) 



(q) Ui(T), (25) 



where 



f ij= / de^(i-o^-(0. (26) 

Jo 

is the eigen functions overlap factor. In Figj2]we plot the first eigen functions evaluated numerically. The Eq. ([25|) 
shows that the input (at the write stage) quantized amplitude of the signal field in the ith eigen mode gets restored 
at z — L with the factor ~ exp (— i<7 2 L/2fc ) A|/jj. Like in the Raman-type schemes, our model demonstrates the 
time reversal symmetry between the input and output modes, as follows from the input (|24p and output (|25[) signal 
decomposition. The read-in coupling parameter k sa 4 and a read-out coupling parameter in excess of 20 is required 
to achieve memory efficiency > 0.95 [181 ]. 

The storage capacity, that is the number of transverse modes that the quantum hologram can store is one of the 
most important parameters for the memory applications. The memory capacity of the previously proposed [10] thin 
quantum hologram is strongly limited by diffraction. Namely, the diffraction spread of an input image element of the 
linear size d should be small, X/d ■ L < d. Hence, the number TV of the resolved by the thin hologram image elements 
(modes) is determined by the sample Fresnel number Fjy, so that N ~ S/d 2 ~ S/XL — Fjf, where S is the sample 
cross section. The authors of [llj who have considered the multiple optical modes storage in an atomic ensemble with 
appropriately phase-matched backward propagating retrieval also estimate the number of stored transverse modes to 
be equal to F N . 

In comparison, the effect of diffraction on the signal wavefront stored in the volume hologram is the same as in free 
space and can be compensated simply by using a lens system with the input focal plane at z = 0. The limitation on 



7 



the capacity for the volume hologram comes first from the paraxial approximation X/d < e, where e <C 1 is the small 
parameter of the paraxial approximation. The second limitation is due to the geometry of the sample, X/d- L ss ^/S, so 
that the quantized waves propagate inside the atomic cell. Hence the number of stored modes in the volume hologram 
is estimated to be S/d 2 ~ min{e 2 S'/A 2 , F^}. For not too elongated samples with VS/L > e the number of modes is 
equal to e 2 S/X 2 which is larger than for memories proposed so far. There are indications that for the co- propagating 
write-in and read-out A- scheme based multimode memories the capacity is also of the order of e 2 S/X 2 [2Cl |. 



III. CONCLUSIONS 



We have presented an extension of the classical volume hologram into the quantum domain. The volume quantum 
hologram can store entangled and other quantum images and offers a new approach to a parallel spatially multimode 
quantum memory for light. The equations of motion for the quantized field and spin coherence waves are derived and 
examined in paraxial approximation. Our analysis could be of use for the spatially-multimode generalization of other 
quantum memory models. 

The counter-propagating geometry of the classical and quantized fields, combined with the spins rotation in a 
constant magnetic field, allows for the efficient quantum image transfer between light and matter in a single pass 
of light. The volume hologram proposed here does not require any extra operations, such as squeezing, in order to 
achieve, in principle, perfect performance. 

In our scheme, different spatial modes of the incoming field are stored in the corresponding orthogonal spatial modes 
of the atomic ensemble. We show that the quantum volume hologram is less sensitive to diffraction and is capable 
of storing more spatial modes as compared to the previously proposed 1 101 ] thin hologram in the co-propagating 
geometry. 

Although we considered spin 1/2 atoms, our analysis can be easily generalized to an arbitrary angular momentum 
states in alkali atoms provided that the optical detuning is greater than the excited state hyperfine structure [l[. 
Future work on volume quantum holograms will include modified geometries of the interacting waves which should 
allow for relaxing of the requirement of motionless atoms, as well as studies of multimode quantum entanglement 
between light and matter. 
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Commission FP7 under the grant agreement n 221906, project HIDEAS. The authors also acknowledge the support 
of the Russian Foundation for Basic Research under the projects 05-02-19646, 08-02-00771, and 08-02-92504 (D.V.), 
and the support by INTAS (Grant 7904). Part of the research was performed within the framework of GDRE "Lasers 
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